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ABSTRACT 

We determine the possible masses and radii of the progenitors of white dwarfs in binaries from 
fits to detailed stellar evolution models and use these to reconstruct the mass-transfer phase 
in which the white dwarf was formed. We confirm the earlier finding that in the first phase 
of mass transfer in the binary evolution leading to a close pair of white dwarfs, the standard 
common-envelope formalism (the a-formalism) equating the energy balance in the system 
(implicitly assuming angular momentum conservation), does not work. An algorithm equating 
the angular momentum balance (implicitly assuming energy conservation) can explain the 
observations. This conclusion is now based on ten observed systems rather than three. With 
the latter algorithm (the 7-algorithm) the separation does not change much for approximately 
equal mass binaries. Assuming constant efficiency in the standard a-formalism and a constant 
value of 7, we investigate the effect of both methods on the change in separation in general 
and conclude that when there is observational evidence for strong shrinkage of the orbit, the 7- 
algorithm also leads to this. We then extend our analysis to all close binaries with at least one 
white dwarf component and reconstruct the mass transfer phases that lead to these binaries. In 
this way we find all possible values of the efficiency of the standard a-formalism and of 7 that 
can explain the observed binaries for different progenitor and companion masses. We find that 
all observations can be explained with a single value of 7, making the 7-algorithm a useful 
tool to predict the outcome of common-envelope evolution. We discuss the consequences of 
our findings for different binary populations in the Galaxy, including massive binaries, for 
which the reconstruction method cannot be used. 
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1 INTRODUCTION 

Now that double white dwarfs are di scovered regularly (e.g. lMarshl 
ll995ll200rilNapiwotzki et all20oil) it has become more and more 
clear that most of them have a mass ratio close to un ity (e.g. 
iMaxted & Marshlll999tlMaxted. Marsh. & Moranll2002ch . This is 
contrary to what is expect ed from sta ndard populati on synth esis 
calculations fe.g. llben. Tutukov. & Yungelsorll997tlHaiJl998h. A 
possible resolution of this issue was investigated by Nelemans et aj. 
fcood l 2001b). In the first paper the observed masses of three dou- 
ble white dwarfs and the well known core-mass - radius relation 
were used to reconstruct the evolution of the binary back to two 
main-sequence stars. It followed that the first phase of mass transfer 
could not be described by the standard common-envelope formal- 
ism (based explicitly on energy balance, assuming angular momen- 
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turn conservation implicitly), nor by stable Roche-lobe overflow. 
Recent calculations using a detailed stellar evolution code have 
confirmed this conclusion (Van der Sluys et al., in preparation). Sta- 
ble Roche-lobe overflow leads to final double white dwarfs with a 
mass ratio larger than one (e.g. llben et alll997llHarll9 98) and the 
observed masses can only be reached by stars with initial masses 
between about 2.3 and 3.5 M© that fill their Roche-lobes within a 
very small initial separation interval (in order to start mass transfer 
in the Hertzsprung gap ). For stan dard population synthesis assump- 
tions fe.g. lNelemans et all20 04) this interval only accounts for 0.3 
per cent of objects forming white dwarfs, so is inconsistent with the 
observation that about 10 per cent of white dwarfs are close pairs 
iMaxted & Marshll 19991) . Instead lNelemans et ail feOOOl) proposed 
an empirical algorithm based explicitly on angular momentum bal- 
ance (implicitly assuming energy conservation) with a single free 
parameter and concluded that all the observed systems could be 
explained with the same value of the free parameter. The second 
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paper showed that using this algorithm a satisfactory model for the 
Galactic population of double white dwarfs can be obtained. 

Since then, quite a few more double white dwar fs have bee n 
discovered. In particular the SPY project l Napi wotzki et all200lh . 
a large survey on the ESO Very Large Telescope, to measure radial 
velocity variations of some thousand white dwarfs in order to de- 
tect duplicity has, and will, enlarge the known double white dwarf 
sample. We therefore repeat the analysis o flNelemans et al|l2000l) . 
including the new discoveries (Section|2j. Furthermore, we study 
the difference between the standard and alternative method in some 
detail (Section[3j. We then extend the analysis to all binaries with 
at least one white dwarf component (Section [4j and sdB binaries 
(Section |5J in order to determine what the free parameter in the 
alternative method must be to explain the observations. We then 
continue with a discussion of the consequences of the our results 
for the different binary populations (Section^ and round off with 
our conclusions. 



2 RECONSTRUCTION OF THE EVOLUTION OF 
DOUBLE WHITE DWARF BINARIES 

We sta rt with a short revision of the method used in Nelema ns et alJ 
(2000). The fact that observed white dwarfs in binaries were the 
cores of the giant stars from which they descend makes it possible 
to reconstruct the properties, in particular the radii, of these giants. 
On the assumption that the observed white dwarf mass is close to 
the mass of the core of the giant at the onset of mass transfer (i.e. 
that the mass-transfer proceeds on a short time-scale compared to 
the nuclear evolution time), the exact evolutionary phase (and thus 
mass and radius) of the giant at that instant can be obtained for each 
possible initial progenitor mass from single star evolution models. 
For a Roche-lobe filling giant its radius, together with its mass and 
that of the companion, determine what the orbital separation at the 
onset of mass transfer. By comparing this with the orbital sepa- 
ration after the mass transfer, the effect of the common-envelope 
phase on the orbit can be reconstructed. 

In iNelemans et alJ 12000]) only double helium white dwarfs 
were considered and a simple core-mass - radius relation for giants 
with degenerate helium cores was used to reconstruct the properties 
of the giants. Here we take a more general approach and use fits to 
detailed stellar evolution calculations to find all the possible giant 
stars that have a core with a mass equal to th e observed w hite dwarf 
mass. For this we use the Hurlev, Pols, & Tout 1 2000) fits which 
enable us to use all observed double white dwarfs, independent of 
them being helium or (low-mass) carbon-oxygen white dwarfs. 

Our exact procedure is as foll ows. For an observed white 
dwarf mass Mwd we use the H urlev et al J 12000) equations to cal- 
culate the masses M g i ant and radii i? g i an t of all the giants which 
have exactly such a core mass. We do this for initial masses of 1, 
1.1, 1.2, ...Mo up to the mass for which the initial core mass, at the 
end of the main sequence, is larger than the observed white dwarf 
mass. While evolving the stars we keep track of the maximum ra- 
dius the star has reached previously so that only giants that actu- 
ally can fill their Roche lobe are selected. Finally, we only consider 
stars if they have passed through the Hertzsprung gap and have de- 
veloped convective envelopes. Radiative stars in the Hertzsprung 
gap can avoid a common-envelope phase so that our assumption of 
mass transfer on a short time-scale compared to the evolutionary 
time-scale is not appropriate. 

For each of the possible masses m for the companion (see be- 
low) we use the size of the Roche lobe Rl in units of the separation 



Table 1. Properties of the observed double white dwarfs 



Object (WD/HE) 


P 


A4\VD,2 


-M\VD,1 


Ref 




(d) 


(M ) 


(M ) 




0135-052 


1.56 


0.47 


0.52 


1,2 


0136+768 


1.41 


0.47 


0.37 


8,13 


0957-666 


0.06 


0.37 


0.32 


3,7 


1022+050 


1.16 


0.35 




8 


1101+364 


0.15 


0.29 


0.35 


4,13 


1115+116 


30.09 


0.52 


0.43 


12 


1202+608 


1.49 


0.4 




6 


1204+450 


1.60 


0.46 


0.52 


8,13 


1241-010 


3.35 


0.31 




5 


1317+453 


4.87 


0.33 




5 


1349+144 


2.12 


0.44 


0.44 


14 


1414-0848 


0.518 


0.71 


0.55 


11,15 


1428+373 


1.143 


0.33 




9 


1704+481 


0.14 


0.39 


0.56 


10 


1713+332 


1.12 


0.35 




5 


1824+040 


6.27 


0.39 




8 


2032+188 


5.084 


0.36 




8 


2209- 1444 


0.28 


0.58 


0.58 


16 


2331+290 


0.17 


0.39 
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Refere nces: (11 [^affe^ Liebert & Olszewski ll988|): (21 iBer geron et alJ 
ll989h: (3) iBragagLta et alJ ll99dt) : (41 iMarshl Il995j): (5) 
Marsh. Dhillon. & DucJd ll995t): (6) iHolberg et alj 1 1993): (7) 
Moran. Marsh. & BragagliJ <1997l) : (8) iMaxted & Marshl <1999l): (91 
M^n T200oFand P. Maxted, pri v ate com munica tion: (101 iMaxtedet alj 
1 2000a ): (111 INapiwotzki et alj <2002ah : (12 1 IMaxted et alj <2002aF : 
(131 IMaxted. Marsh. & Moran) l2002d): (14) iKarl et alj 1200211 : (151 
INapiwotzki et alj<2002bl) : (16l lKarl et all <2003l) : 



a, tl = Rl /a,, as given by Eggleton ( 1983) to determine the sep- 
aration at the onset of the mass transfer assuming i? g i an t = Rl- 
The range of companion masses considered is determined by the 
observations. If the mass of the companion is known that mass is 
used but for unseen companions in double white dwar f systems we 
use the extremes of 0.2 and 1.4 Mq as in lNelemans et alj j2000). 

In Table □ we list the properties of the observed double 
white dwarfs. It includes both updates and additions to table 1 of 
INelemans et alj l2000l) . There are now 10 binaries in which the 
masses of both components are known. For these we can use our 
reconstruction method twice, first for the last phase of mass trans- 
fer in which the white dwarf with mass Mwd, 2 is formed and the 
companion star was a white dwarf of mass Mwd,i- This gives the 
separation before the second phase of mass transfer and the mass of 
the giant that formed white dwarf 2. We then calculate the separa- 
tion after the first phase of mass transfer by assuming the separation 
only changed owing to mass loss in a wind from the progenitor of 
white dwarf 2. Finally we use the initial mass of the progenitor of 
white dwarf 2 and the mass of white dwarf 1 (Mwd,i) to calcu- 
late the change in separation in the first phase of mass transfer. The 
only extra constraint we have to put in is that we require the pro- 
genitor of white dwarf 1 to be more massive than the reconstructed 
progenitor of white dwarf 2. 

We now discuss the results for the first phase of mass trans- 
fer because that is the phase that was found to be incons istent with 
the sta ndard common-envelope formalism, proposed by Paczvhski 
( 1976) to explain the existence of short-period binaries with white 
dwarf components and cataclysmic variables. It is generally as- 
sumed that the outcome of the common-envelope phase is deter- 
mined by the energy balance, implicitly assuming angular momen- 
tum conservation. I.e. that the orbital energy of the binary is used 
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Figure 1. Left: Reconstructed 7 values for the first phase of mass transfer in the formation of double white dwarfs. Right: reconstructed aX values for the 
same. The horizontal lines are made up of small dashes representing the reconstructed values of 7 and aX for different values of the mass of the progenitor 
of the white dwarf and the companion. The different lines for each object represent different values of the white dwarf mass (within 0.05 Mq of the value in 
Tabled. 



to expel the en velope of the giant with some efficiency a (e.g. 
Webb inkll98l 



GM & Me (GM c m 



2a f 



GM s m 
2ai 



(1) 



where subscripts g, e and c are for giant, envelope and core 
respectively and we assume the companion mass does not 
change during the common-envelope phase. The structural pa- 
rameter A is normally taken as a con stant ( e.g. A = 0.5 
Ide Kool. van den Heuvel. & PvlvseJll987h . or as iNelemans et alJ 
(2000) the A factor is incorporated in the uncertain efficiency factor 
to give one free parameter aA and this is what we do here. We will 
refer to this method as the standard a-formalism. 

The algorithm based explicitly on the equation for angular 
momentum balance (implicitly assuming energy conservation, but 
not necessarily only f or orbital and binding energy) proposed by 
INelemans et alJ OOOO) is described by 



AJ 
J 



AM t , 



1 



M s + m' 



(2) 



In the remainder of the paper we refer to this method as the 7- 
algorithm. 

For each double white dwarf for which the masses of both 
white dwarfs are known we can calculate the range of possible 
masses of the secondary from Mwd,2 and separations after the 
first phase of mass transfer and the possible masses and radii of 
the primary from Mwd.i and thence the separation at the onset of 
the first phase of mass transfer. That means that all terms in equa- 
tion except aA and all terms in equation J2j except 7 are known. 
For the calculation of the total angular momentum we include the 
angular momentum of the giant, assuming all the angular momen- 
tum resides in the envelope, which we approximate as an n — 3/2 
polytrope. In Fig.Qwe show, for each of the observed double white 
dwarfs, the possible values of aA and 7 that we find in this way. 
Each possible combination of progenitor and companion mass is 
shown as a small dash, forming horizontal lines. The different lines 
for each object are for different values of the white dwarf mass 
to account for measurement errors (which we take as ±0.05 Mq). 
WD1 1 15+1 16 is shown twice because it is not clear from the obser- 
vations which of the two white dwarfs is white dwarf 1 and which 
is white dwarf 2. 




M giant 

Figure 2. Reconstructed 7 values versus the mass of the giant for the first 
phase of mass transfer in the formation of double white dwarfs. 



We confirm the findings o f Nele mans et alJ feOOOh that the first 
phase of mass transfer in the evolution leading to the observed dou- 
ble white dwarfs cannot generally be described by the standard 
a-formalism because the reconstructed values of aA are negative. 
The only exception is WD1704+481 which does have a mass ra- 
tio in the range expected from evolution governed by the standard 
a-formalism. 

As to the value s of 7, we also recover the results of 
INelemans et alJ J2000h with typical values around 1.5 with a large 
spread. However for a value of 7 between 1.5 and 1.75 we can find 
simultaneous solutions for all objects. 

To asses the likelihood of the solutions found with 7 around 
1.5 we plot the mass of the giant versus the reconstructed 7 value 
in Fig.fi] Typical giant masses are between 1.5 and 2 Mq, just as 
one would expected for the more massive components in binaries 
that eventually form double white dwarfs, i.e. in which both stars 
evolve off the main sequence within the age of the Galaxy. 
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Figure 3. Relative change in separation for dynamical mass transfer described by the standard a-formalism (left panel) and the 7-algorithm (right panel) as 
function of mass ratio q and core mass fraction /i. The logarithm of (a;/af ) is shown as the grey scale and the contours. 



3 COMPARISON OF THE STANDARD a-FORMALISM 
AND THE 7-ALGORITHM 

lNelemanset"ail l2000) proposed that the first phase of mass trans- 
fer in the evolution to a close double white dwarf was special in the 
sense that it is most likely a phase of dynamical mass transfer but 
in a binary with mass ratio not too far from unity. In such a binary 
the angular momentum of the orbit is so large that the envelope of 
the giant can be spun up easily. This removes the drag forces that 
might drive any loss of orbital energy. We will discuss the question 
of the physical interpretation and applicability of the 7-algorithm in 
a forthcoming paper but here we consider the effect of both meth- 
ods described above on the change in orbital separation for a wide 
range of giant and core masses and mass ratios. 

The change in separation of the binary for the standard a- 
formalism is 

(-) =^( i+ ^r> (3) 

Vai/ a M g \ ctXr^mJ 

where we used R g = i?L, giant = r^ai and we have again assumed 
that none of the envelope is accreted by the companion. For the 
7-algorithm this ratio is 

/OA _ (M g m\ 2 ( M c +m \( M e V 

\aJ J \M c m) \M s +mJ\ 1 M g + rnJ ' U 

Only ratios of the masses of the different components (gi- 
ant, companion, core, and envelope) enter these equations, so that 
the relative change in the orbital separation does not depend on 
the total mass in the system but only on the mass ratios of the 
different components. There are only three independent masses 
(companion mass, giant mass and either core mass or giant en- 
velope mass). These are characterised by only two ratios and the 
ratios q — M s /m and p = M COTC /M g or alternatively A = 
M c /A'fgiant = 1 — H conveniently simplify equations 13141 : 




With these equations we calculate the change in separation for 



both methods as functions of q and fi, using aX = 2 and 7 = 
1.5. The results are shown in Fig. [3] The grey scale denotes the 
log of af/at, with darker shades a greater shrinkage of the orbit. 
Contours of constant logaf/ai of 0.5, 0, —0.5 etc are also shown 
in the figures. Fig.|3|shows that for any combination of q and /i the 
standard a-formalism gives a strong shrinkage of the orbit while 
with the 7-algorithm there is a wide range from expansion to very 
extreme shrinkage (or even guaranteed merger if all the angular 
momentum is lost). 

The reconstruction of double white dwarfs discussed in Sec- 
tion[2|finds evidence for a strong reduction in the last but not in the 
first phase of mass transfer. For typical progenitors of double white 
dwarfs, with n between 0.2 and 0.5 and q between 1 and 2 in the 
first and between 2 and 6 for the second phase of mass transfer the 
standard a-formalism would give strong shrinkage of the orbit in 
both cases. In contrast, the 7-algorithm gives widening or very mild 
shrinkage in the first, but strong shrinkage in the second phase of 
mass transfer and thus might explain both phases in the evolution 
to double white dwarfs. 

We explore the difference between the two methods further in 
Fig.|4]where, for initial masses of 1, 2, and 3 M©, we plot a grey 
scale of the period at the end of the common-envelope phase for 
both methods as function of the core mass and companion mass. We 
again use 7 = 1.5 and aX = 2. It can be seen that for the standard 
a-formalism the final periods are below 10 d except for the most 
massive cores, while for the 7-algorithm it depends strongly on the 
mass of the companion. Indeed, for relatively high core masses and 
companion masses, very large final periods, above 1000 d, can be 
induced. This is interesting in the light of the existence of symbi- 
otic binaries, barium and S-stars with periods in that range. The 
alternative for the formation of these binaries is that they avoided 
a common-envelope phase. That is they have stable mass transfer 
or avoid mass transfer at all and would be expected to have even 
longer orbital periods. We will come back to these binaries in Sec- 
tion|7j 

The observational requirement for strong orbital shrinkage has 
always been for rather extreme-mass-ratio systems such as cata- 
clysmic variables and low-mass X-ray binaries. So the fact that the 
7-algorithm actually produces a strong shrinkage at large mass ra- 
tios makes it useful to consider the 7-algorithm in more extreme 
mass ratio common-envelope phases, such as the last phase of mass 
transfer leading to a close double white dwarf and mass transfer in 
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Figure 4. Final periods as function of core mass and companion mass with the 7-algorithm (left) and the standard a-formalism (right) for giants of 1 Mq 
(top), 2 Mg (middle) and 3 Mq (bottom). The darkest shades represent periods of 0.01 d and the lightest periods of 1000 d. Dashed contours are for constant 
log (P/d) = -2, -1.5 -0.5 and the solid contours are for log (P/d) = 0, 0.5, 2.5. The even light gray area shows parts of parameter space for which stable 
mass transfer is expected. The gap in the middle occurs because the core mass grows during stages (core helium burning) when the star has a smaller radius 
than it had before and Roche-lobe overflow cannot take place. The white area below the shaded areas denotes combinations for which the systems merge when 
using the 7-algorithm because all the angular momentum is lost from the system. 



binaries leading to a close binary with a white dwarf and a main- 
sequence star, many of which are observed. We can use the same 
procedure we used to reconstruct the first phase of mass transfer in 
double white dwarfs to reconstruct any of these. 



4 WHITE DWARF BINARIES 

There are two classes of objects for which the standard a- 
formalism has been successfully used to explain their properties. 
These are the last phase of mass transfer leading to the forma- 
tion of a close double white dwarf (e.g. iNelemans et al] f2001b) 
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Figure 5. Left: Reconstructed 7 values for the last phase of mass transfer in the formation of double white dwarfs (see Table HI. Right: reconstructed a A 
values for the same. 



and the formation of close white dwarf - main-sequence binaries. 
The latter ar e expected to be the precursors of cataclysmic variables 
(e.g. de Kool & Ritter 1993). We discuss these binaries here in turn, 
comparing again the standard a-formalism and the 7-algorithm in 
order to asses how well they do in predicting the outcome of the 
common-envelope phase. 

There are two complicating factors which we have to take into 
account, the first of which is tidal interaction. If the spin angular 
momentum of one of the components in a binary exceeds one third 
of t he orbital angular momentum, the tidal interaction is unstable 
(see lHulll980llSparks & StecheiH 19741) showed that for mass ra- 
tio's larger than about 6, the tidal instability sets in before the giant 
fills its Roche lobe. Some of the binaries we shall discuss must 
have had quite extreme mass ratios at the onset of the mass transfer 
because the companions are either low-mass white dwarfs or low- 
mass main-sequence stars. We therefore build in a check for tidal 
stability in the reconstruction process. When a progenitor system 
is found to be tidally unstable for our assumption that the onset of 
the mass transfer is caused by Roche-lobe overflow, we relax this 
assumption and instead assume the mass transfer was caused by 
the tidal instability and we calculate the initial separation at which 
the instability sets in at exactly the right core mass. The last aspect 
of the new procedure is a check whether the initial separation is 
small enough that the companion will actually keep the giant in co- 
rotation with the orbit, because otherwise the tidal instability will 
not set in at all and ma ss transfer is avoided. We u se the max imum 
separation as given in lNelemans & Taurisl Il998t) . based on IZahnl 

The second is the question whether the current orbital period 
of the observed systems is a good estimate of the post-mass-transfer 
period. In particular the systems with a low-mass main-sequence 



companion might have experienced angular momentum loss owing 
to magnetic brak ing 1 Verbunt & Zwaaall9811) . In a recent study 
ISchreiber & Gan sicke (2003) carefully investigated this effect in 
30 post-common-envelope binaries and found that virtually all ob- 
served periods are close to the initial periods after the common en- 
velope. Only for EC 13471-1258 and BPM 71214 did they find any 
evidence for significant orbital evolution. Even in these cases the 
change is relatively small, so for the current purpose we prefer to 
use the observed periods in the analysis. 

4.1 Double white dwarf binaries 

For double white dwarfs we can reconstruct the last phase of 
mass transfer for all 19 objects listed in Table Q We find that 
the reconstructe d values of a\ are indeed in a reasonable range 
(as in lNelemans et alj2 000). Most systems can be explained with 
qA ~ 0.5. However the spread is large. As for the reconstructed 
values of 7 we again find that all systems can be explained with a 
value of 7 ps 1.5. All reconstructed values of 7 and ot\ are shown 
in Fig. [5] As before WD1 115+1 16 is included twice because it is 
unclear which of the two objects is formed last. 

4.2 Pre-cataclysmic variables and other white dwarf - main 
sequence stars 

For the properties of the observed pre-cataclysmic variables and 
other white dwarf - m ain-seque nce binaries w e use the compila- 
tion of iHillwig. Honevcutt. & Robertson! ^2000) extended and up- 
dated with recent published results and systems not in their table. 
All details are given in the appendix, in Table lAll Most objects are 
short-period systems in which the companion to the white dwarf 
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Figure 6. Left: Reconstructed 7 values for the mass transfer in white dwarf, M dwarf (or earlier type) star binaries (see Table I A lV Right: reconstructed a\ 
values for the same. 



is a low-mass main-sequence star. Towards the bottom of the ta- 
ble (which is sorted by increasing orbital period) there are a few 
interesting systems that have rather large orbital periods and so 
w ould b e difficult to explain with the standard a-formalism. Indeed 
in lNelemans et alj feOOd) S1040 and AY Cet were cited as further 
evidence for the 7-algorithm. 

The results of the reconstruction give a quite similar pat- 
tern to the last phase of mass transfer in the evolution of double 
white dwarf binaries: both methods can more or less explain all 
the observed systems. All values are shown in Fig. [6] As men- 
tioned above, at the long-period end (top of Fig. |6|, the standard 
Q-formalism cannot explain a few systems. The values of aX also 
seem to correlate with the final periods: the lower half of Fig. [6] 
i.e. the shorter orbital periods, requires lower values of aX than the 
upper half. 



5 SUB-DWARF B BINARIES 

The last group of binaries we consider in some detail are binaries 
in which one component is a sub-dwarf B (sdB) star. These are 
thought to be helium-burning stars with a very thin hydrogen en- 
velope jHebeJ fT986y When the core helium burning ceases, they 
are expected to settle on the white dwarf cooling branch in the 
HR diagram. Almost all of them are members of a binary system. 
There are essentially two scenarios to form an sdB star in a binary, 
a giant with a non-degenerate helium core loses its envelope to a 
companion (e.g. lHan et al.l2 002) or a a giant with a degenerate he- 
lium core loses its envelope to a companion just before it reaches 
the core mass at which the helium in the degenerate core ignites 
(proposed first bv lD'Cruz et aljl99d assuming the en velope was 
lost by a strong stellar wind). Detailed calculations (e.g. lHan et alJ 



2002) show that there is a small range of core masses for which the 
latter occurs. The fact that on the first giant branch the radius of 
the giant increases very rapidly with the growth of the core mass, 
means that this small range of core masses corresponds to a large 
range in radii so no fine tuning is needed to get Roche-lobe overflow 
for these core masses. The initial reasoning for this scenario came 
from the fact that the observed sdB stars seemed all to be incred- 
ible similar in mass, around 0.5 Mq. This would follow naturally 
from the mass at which helium ignites (about 0.47 Mq). Note that 
Han et al. 1 2002) showed that the mass at helium ignition for stars 
initially above 1.5 Mq drops quite strongly to about 0.33 Mq at 
stars of initial mass 2 Mq as the flash becomes less and less degen- 
erate, but this depends critically on the assumed core overshooting. 

Because sdB stars are bright, they are relatively easy to study 
and, in the las t few years, a large fraction has been surveyed for 
dupli c ity (e.g. iKoen et alJll998t i Moran et alJll999t iMaxted et alJ 
l200lt Morales-Rueda et al. 2003). Many turn out to be close bi- 
naries and for some the mass of the companion can also be deter- 
mined. The properties of these binaries are given in the appendix in 
Table lA2l We note that the general assumption is that all sdB stars 
have a mass of 0.5 Mq, so only entries with a different mass have 
actually been determined in detail. In Fig. we show the recon- 
structed values of aX and 7. We find again that most systems can 
be explained either by a rather low value of aX or a value of 7 close 
to 1.5 The one wide binary (HD 185510 with an orbital period of 
20.7 d) only has a solution for the 7-algorithm. 

The sdB binaries for which the mass of the companion is not 
known are listed in Table lA3l We only list the period and the min- 
imum companion mass derived from the period and velocity am- 
plitude of the sdB stars with the assumption that the binary is seen 
edge on. For the sdB stars with no independent determination of 
the mass we use the results of lHan et alJ y002) and consider all 
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Figure 7. Left: Reconstructed 7 values for the mass transfer in the sdB binaries in which the mass of the companion is known (see Table lA2l . Right: 
reconstructed aX values for the same. 



core masses that are just about to ignite helium plus all helium core 
masses that are non-degenerate as possible progenitors of the sdB 
stars. Since this implies such a wide range of progenitors, almost 
all values of aX and 7 are possible. In Fig. lAll we show the recon- 
structed values. We use as limits to the compan ion mass the min- 
imum mass, derived from the mass function I Morales-Rue da et alJ 
l2003h and 1.4 Mq because most are expected to be white dwarfs. 
If they are main- sequence stars their mass should be much lower. 



6 STATISTICS 

Before considering the physics of the common-envelope process, 
we have a look at the statistics of the reconstructed 7 and aX val- 
ues shown in the figures in the previous sections. Fig. [8] contains 
histograms of the reconstructed values of 7 and aX for the binaries 
we have considered. It is normalised in such a way that the sum of 
all bins of the histograms is equal to the total number of systems we 
reconstructed. The distribution of 7 values is actually peaked at 1.5 
while the aX distribution is skewed to lower values. This suggests 
that the 7-algorithm is a useful tool for predicting the outcome of 
common-envelope evolution. Another surprising result is that each 
of the distributions looks very similar for the different types of bi- 
naries. This suggests that the processes determining the mass trans- 
fer outcome are similar for the different types and so the different 
types of binaries can be seen as independent measurements of the 
same process. Furthermore it shows that the outcome of the com- 
mon envelope apparently does not depend strongly on whether the 
companion star is a very compact white dwarf or a much larger 
main-sequence star. 

Some caution is needed in simply comparing the numerical 
values of aX and 7. The definition of 7 implies a certain restriction 
to the values that can be obtained by the reconstruction method: 
AJ/J has a value between and 1, and M e /(M S + m) can be 
written as M c /(M g (l + <?)), with M c /M s limited between 0.4 and 
0.9. Combining these numbers, we simply never can find values of 

7 outside the range to 5. A slightly smaller range is expected for 
very short period systems, where the final angular momentum is 
small. That means that AJ/J is or the order unity, leading to val- 
ues of 7 larger than 1 . On the other hand, the standard a-formalism 
is less constrained. Owing to its definition (-Ebind/ AE^) a can 
never be zero because the binding energy never is zero. Another 
property of the standard a-formalism is that for most cases A£" rb 



is completely determined by the final separation. The result of that 
is that if the final separation is not determined by the binding en- 
ergy, giants with similar binding energy will show a reconstructed 
a which is correlated with the final period, which seems to be the 
case in Fig.|6| 

We have deliberately kept the discussion focused on the obser- 
vations and tried to interpret them as model independently as possi- 
ble. The only theoretical ingredients so far are evolution models of 
single stars. Though the fact that all binaries can be explained with 
a single value of 7 gives a useful tool to describe the outcome of a 
common-envelope phase, it doesn't give a physical understanding 
of the process. Our ideas about the interpretation of the 7-algorithm 
in terms of a possible physical mechanism will be discussed in a 
forthcoming paper (Nelemans & Tout, in prep.). 



7 CONSEQUENCES OF THE 7-ALGORITHM 

Although we showed above that all observed double white dwarfs, 
pre-CVs and sdB binaries could be explained with the 7-algorithm, 
there are some drastic consequences of this assumption. The most 
important is that the simple form of equation immediately 
shows that for given 7 there is a limit to the amount of mass that 
can be lost before the system merges. E.g. for 7 = 1.5, the system 
will lose all its angular momentum for M e > 2/3(A/ g + m). For 
extreme mass ratios this is similar to a core mass fraction /1 < 1 /3, 
which is often realised in the early evolution of stars. For binaries 
which undergo a first phase of stable mass transfer, the secondary 
often accretes enough mass that in the second phase of mass trans- 
fer the system will merge due to the extreme mass ratio. Another 
situation where a large fraction of the total mass is lost is in the case 
of mass trans fer from a giant to another giant leading to a double 
spiral-in (e.g. Brown 1995; Nelemans et al. 2001b). If in that situ- 
ation the 7-algorithm is used, it normally leads to complete merger 
of the two cores. 



7.1 Double white dwarfs and (pre-)CVs. 

In order to assess the consequences of using the 7-algorithm 
in all phases of dynamically unstable mass transfer, we made a 
population synthesis cal culation with a ll assumptions ident ical to 
the model described in Nelema ns. Yungelson. & Portegies Zwartl 
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Figure 8. Histograms of reconstructed values of 7 and a A for the last phase of mass transfer in the formation of double white dwarfs (top left), in white dwarf 
plus M dwarf binaries (bottom left) and the two kinds of sdB binaries (top end bottom right). 



(2004), except for the common-envelope phase, where we now use 
the 7-algorithm in all cases. 

For double white dwarfs the main effect is that the number of 
observable systems in the Galaxy goes down by almost a factor of 
3. On the one hand this is due to the fact that more systems merge 
because all the angular momentum is lost. This happens for systems 
in which the secondary is rather massive by itself or because it ac- 
cretes in the first phase of mass transfer, or because mass transfer 
begins when both stars are giants. On the other hand more systems 
form with such long orbital periods (above 40 d) so that current 
observing programs are not sensitive to them. This happens for the 
the systems forming from low-mass stars. 

The period distribution of the model in which we use the stan- 
dard a-formalism in the last p hase of m ass transfer matches the ob- 
served distribution quite well (N elemans et alj2001bl) so we would 
expect the new model to do less well. However, this is not signifi- 
cantly the case. The reason is that, because of the relatively small 
number of observed systems, we are essentially comparing the ob- 
served period range with that in the model. The short end of this 
range is largely determined by the fact that systems merge due to 
angular momentum loss by gravitational wave radiation and thus 
disappear from the observable sample and the long period end by 
the limits of the current methods of period determinations. 



The difference in the model is very important for the merger 
rate of double white dwarfs. It drops by about a factor 6. For mas- 
sive pairs the situation is even more dramatic. The merger rate 
of pairs with a mass above the Chandrasekhar mass reduces from 
1.1 x 10~ 3 yr _1 to 7.7 x 10 _6 yr _1 . Similarly, the confusion lim- 
ited noise background of the Galactic double white dwarf popula- 
tion for the space based gravitational wave detector LISA decreases 
by about a factor 2. 

An even greater change is that the mass-ratio distribution of 
the close double white dwarfs is even more peaked around unity 
when the 7-model is used throughout. This means the chances 
for double white dwarfs to start stable mass transfer and evolve 
into AM CVn systems are significantly reduced. The birth rate of 
AM CVn systems drops by a factor 20 from 1.3 x 10~ 3 yr _1 to 
6.4 x 10 -5 yr -1 , while the total number of systems in the Galaxy 
goes from 2.3 x 10 7 t o 9.2 x 10 5 . For AM C Vn systems formed 
from helium stars (see lNelemans et alll2001a|. for a discussion of 
the ways to form AM CVn systems) the reduction is even larger, 
about a factor of 100. With such a small number it would become 
problematic to explain the number of known systems which are be- 
lieved to be only a small part of the total observable population (e.g. 
iNelemans et all2001 si . 

As expected from the fact that we can explain most observed 
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white dwarf - main-sequence binaries, this population is not much 
affected by the use of the 7-algorithm. Of course the white dwarfs 
with relatively massive companions have longer periods compared 
to the case when we use the q-for malism in the first phase of mass 
transfer (see lNelemans et all200fj) . There is quite a large effect on 
the formation of cataclysmic variables, with their current birth rate 
reduced by a factor 2.5, but the total number in the Galaxy by a fac- 
tor 7 (from 26.8 to 3.8 million). However, these numbers are within 
the range expected from observations (e.g. Marsh 2001), particu- 
larly because, in our Galactic model, the fraction of systems close 
to the sun is lower than in an exponential disc (see Ne lemans et alJ 
2004). 

7.2 Consequences for other binaries 

Some of the properties of the 7-algorithm are relevan t to the for- 
matio n of symbiotic stars, barium and S-stars (e.g. ljorissen et all 
Il998h . As shown in Fig.|4| it can lead to quite wide binaries after 
unstable mass transfer, otherwise the observed long periods of bar- 
ium stars mean that population models rely on either wind accre- 
tion t o transfer s-process enriched materi al (e.g. lBoffin & Jorissenl 
Il988t iKarakas. Tout. & Lattanziol 12000) or extra mass loss on 
the AGB to avoid a common envelope jTout&EggletorJll988t 
lHan et all 19 95). A detailed anal ysis of the post- AGB binary in the 
Red Rectangle iMen'shchikov et al. 2002) suggests an evolution- 
ary scenario in which the 7-mechanism is needed to explain the 
current system parameters. 

Similar problems also affect higher-mass binaries, particularly 
the common-envelope phases in the evolution leading to low-mass 
X-ray binaries and double neutron stars. The standard scenario for 
the formation of low-mass X-ray binaries Ivan den Heuvelll983h 
involves the common-envelope evolution of a star that, after los- 
ing its envelope, becomes a neutron star or even black hole and a 
low-mass main-sequence star. A relatively low-mass neutron star 
progenitor with an initial mass of 9 M©, which attains a maximum 
core mass of abo ut 2.5 Mp, w hen its total mass is 8.5 M©, accord- 
ing to the|Hurie^^|200(3) formulae, has M c /(M g +m) > 0.63 
and so must have 7 < 1.58 in order not to merge. It turns out that 
for more massive stars this limit on 7 decreases only rather slowly, 
so for values of 7 not too much in excess of 1.5 formation of low- 
mass X-ray binaries is still possible. For the formation of double 
neutron stars the situation is quite similar. 



8 CONCLUSIONS 

We have used the masses of observed white dwarfs in binaries to 
estimate the radii of their progenitors on the assumption that the 
white dwarf masses are good approximations to the core masses of 
their progenitor giants. Using these progenitor masses we have re- 
constructed the parameters of the progenitor binary systems of the 
observed white dwarf binaries. These, together with the observed 
binary parameters, were used to reconstruct the change in orbital 
separation during the mass transfer phase in which the white dwarf 
was formed. By comparing this change to the expected change for 
the standard Q-formalism, explicitly based on the energy balance, 
and an the 7-algorithm, explicitly based on the angular momen- 
tum balance, we derived the values of the free parameters in these 
methods. 

The main result is that, as was found earlier, for the first 
phase of mass transfer in the evolution leading to the currently 
observed double white dwarf systems, the standard a-formalism 



cannot explain the observations, while the 7-algorithm can. For all 
the other reconstructed phases either method can explain the ob- 
servations. However, the reconstructed values for the 7-algorithm 
strongly cluster around 1.5, while the values of the free parame- 
ter in the standard Q-formalism (the efficiency parameter), show a 
wide range of values, skewed towards low (< 0.5) values. Thus 
the predictive power, at least in statistical sense, of the 7-algorithm 
seems to be greater than the standard Q-formalism. 
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Table Al. Properties of the observed white dwarf, other star binaries 

Object P Af WD m Ref 

(d) (M ) (M ) 

GD448 0.103042 0.41 0.096 1 

NNSer 0.130080 0.570 0.120 1 

EC 13471-1258 0.15074 0.78 0.43 10 

MS Peg 0.173660 0.480 0.220 1 

BPM 71214 0.20162 0.77 0.4 9 

PG 1224+309 0.25869 0.45 0.28 1 

CCCet 0.286654 0.400 0.180 1 

RRCae 0.304 0.467 0.095 1 

BPM 6502 0.33678 0.5 0.15 8 

GKVir 0.344331 0.510 0.100 1 

KVVel 0.357113 0.630 0.25 11 

UUSge 0.465069 0.630 0.290 1 

G1781A 0.497 0.35 0.25 4 

V471 Tau 0.521183 0.84 0.93 2 

HZ 9 0.564330 0.510 0.280 1 

UXCVn 0.573703 0.390 0.420 1 

PG 1026 0014 0.597257 0.65 0.220 15 

EGUMa 0.667650 0.64 0.42 13 

REJ2013+400 0.7056 0.56 0.18 12 

REJ1016-0520 0.789 0.61 0.15 12 

ESTCMa 1.262450 0.57 0.39 12 

BEUMa 2.291166 0.7 0.36 1 

HD33959C 2.99 0.6 1.5 14 

Feige24 4.231600 0.49 0.37 13 

G203-047ab 14.7136 0.6 0.33 3 

V651Mon 15.991000 0.400 1.800 1 

IK Peg 21.721700 1.100 1.700 1 

S1040 42.8 0.22 1.5 6 

AYCet 56^8 025 2^2 7_ 

References: 

(11 IHillwig Honevcutt. & RobertsorJ <200(j>: (2) lO'Brien Bond. & Siorl 
20011) : (31 iDelfosse et alj Il999l): (4) iGiziJ il998l): (5) Benedict et all 
2000[) : (61 Landsman et alj <199tI) : (7 1 ISimon. Fekel & Gibson jj 
19851): (8) iKawkaetalJ COOCl): (91 iKawka et alj 12002^ (10) 
O'Donoghue et alj 120031): (111 Iffilditch. Harries. & Hilll Il99d>: (12) 
Vennes. Thorstensen. & Polomski TT9991) : (131 iBleach et alj fcOOCj) (141 
Vennes. Christian. & Thorstensenl Il998j). note tha t IHill wig et alj <20001) 
list the parameters of HD 33959A; (15l lSaffer etal]|l993h 



Table A3. Periods and minimum companion masses of sdB binaries, from 
Moral es-Rueda et alj <2003t) 



Object 


P 


Mnin 




(d) 


(M ) 


PG1017-086 


0.0729939 


0.066 


PG1043+760 


0.1201506 


0.106 


PG1432+159 


0.22489 


0.294 


PG2345+318 


0.2409458 


0.379 


PG1329+159 


0.249699 


0.083 


TWCrv 


0.328 


- 


PG1 101+249 


0.35386 


0.424 


KPD 1946+4340 


0.403739 


0.628 


PG1743+477 


0.515561 


0.438 


PG000 1+275 


0.528 


0.293 


PG0101+039 


0.569908 


0.370 


PG1725+252 


0.601507 


0.381 


PG1247+554 


0.602740 


0.090 


PG1248+164 


0.73232 


0.207 


PG0849+319 


0.7451 


0.228 


PG1627+017 


0.829226 


0.273 


PG1 116+301 


0.85621 


0.356 


PG0918+029 


0.87679 


0.313 


HE1047-0436 


1.213253 


0.458 


PG0133+114 


1.2382 


0.388 


PG15 12+244 


1.26978 


0.458 


UV01735+22 


1.278 


0.539 


HD 171858 


1.529 


0.510 


PG17 16+426 


1.77732 


0.366 


PG1300+279 


2.2593 


0.346 


PG1538+269 


2.501 


0.600 


KPD0025+5402 


3.571 


0.235 


PG0839+399 


5.622 


0.226 


PG0907+123 


6.1163 


0.521 


PG1032+406 


6.779 


0.247 


PG0940+068 


8.330 


0.634 


PG1 110+294 


9.415 


0.633 


PG1619+522 


15.357 


0.376 


PG0850+170 


27.81 


0.466 



GD448 


0.103042 


0.41 


0.096 


NNSer 


0.130080 


0.570 


0.120 


EC 13471-1258 


0.15074 


0.78 


0.43 


MS Peg 


0.173660 


0.480 


0.220 


BPM 71214 


0.20162 


0.77 


0.4 


PG 1224+309 


0.25869 


0.45 


0.28 


CC Cet 


0.286654 


0.400 


0.180 


RRCae 


0.304 


0.467 


0.095 


BPM 6502 


0.33678 


0.5 


0.15 


GKVir 


0.344331 


0.510 


0.100 


KVVel 


0.357113 


0.630 


0.25 


UUSge 


0.465069 


0.630 


0.290 


G1781A 


0.497 


0.35 


0.25 


V471 Tau 


0.521183 


0.84 


0.93 


HZ 9 


0.564330 


0.510 


0.280 


UX CVn 


0.573703 


0.390 


0.420 


PG 1026 0014 


0.597257 


0.65 


0.220 


EGUMa 


0.667650 


0.64 


0.42 


REJ2013+400 


0.7056 


0.56 


0.18 


REJ1016-0520 


0.789 


0.61 


0.15 


ESTCMa 


1.262450 


0.57 


0.39 


BEUMa 


2.291166 


0.7 


0.36 


HD 33959C 


2.99 


0.6 


1.5 


Feige 24 


4.231600 


0.49 


0.37 


G203-047ab 


14.7136 


0.6 


0.33 


V651 Mon 


15.991000 


0.400 


1.800 


IK Peg 


21.721700 


1.100 


1.700 


SI 040 


42.8 


0.22 


1.5 


AYCet 


56.8 


0.25 


2.2 



Table A2. Properties of the observed sdB binaries, for which the mass of 
the companion is known. 



Object 


P 


m 


M 


Ref 




(d) 


(Mq) 


(Mq) 




PG 1017-086 


0.073 


0.5 


0.078 


9 


KPD 0422+5421 


0.090180 


0.51 


0.526 


4 


KPD 1930+2752 


0.095111 


0.5 


0.97 


7 


HS0 705 +6700 


0.095647 


0.48 


0.134 


8 


PG 1336-018 


0.101 


0.5 


0.15 


3 


MTSer 


0.113227 


0.6 


0.2 


1 


HWVir 


0.116720 


0.48 


0.14 


5 


AA Dor 


0.261540 


0.330 


0.066 


6 


V477 Lyr 


0.471729 


0.51 


0.15 


1 


FF Aqr 


9.207755 


0.5 


2.0 


1 


HD 185510 


20.7 


0.304 


2.27 


2 


References: 



1998); (4) 
20001): (71 



(9) fMaxted et al 



: ( 2lljefferv & Simonl 
Orosz & Wade! ll999t); (5l lWood 



Maxted, Marsh, & North 120001 



. 1 2002b) 
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Figure Al. Left: Reconstructed 7 values for the mass transfer in sdB binaries in which the mass of the companion is unknown. From the radial velocities of 
the sdB star and an assumed sdB mass of 0.5 Mq a minimum mass is inferred. For the upper limit a mass three times the minimum mass is assumed. Right: 
reconstructed qA values for the same. 



© 2004 RAS, MNRAS OOOHTfTTI 



